An instability of the Reissner- Nordstrom solution and 
new hairy black holes in d = b dimensions 

Yves Brihaye,^ Eugen Radu-'- and D. H. Tchrakian** 

^Physique-Mathematique, Universite de Mons, Mons, Belgium 
"'"Institut fiir Physik, Universitat Oldenburg, Postfach 2503 D-26111 Oldenburg, Germany 
, ^Department of Computer Science, National University of Ireland Maynooth, Maynooth, Ireland 

04 I ^School of Theoretical Physics - DIAS, 10 Burlington Road, Dublin 4, Ireland 

> ' 

O ■ November 15, 2011 

^ Abstract 

The d — 5 Reissner-Nordstrom black hole becomes unstable when considered as a solution of Einstein- 

^ ' Yang-Mills-Chern-Simons theory. The existence of a marginally stable mode gives rise to a new branch of 

black holes with non-Abelian magnetic fields outside the horizon. These solutions carry a nonzero electric 
,1 ^ . charge and have finite mass. We argue that these features manifest themselves both for Minkowski and 

j^jq' (Anti-)de Sitter asymptotics. The properties of solutions in a de Sitter background are new to the present 

work and are emphasised. All solutions constructed have finite mass by virtue of the presence of the 
Chern-Simons term, which in d = 5 plays the role of a higher order term scaling appropriately for this 
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purpose. 



1 Introduction 

The Reissner-Nordstrom (RN) black hole is perhaps the simplest and the best known non-vacuum configu- 
ration in general relativity. This solution of the Einstein-Maxwell equations describes the gravitational field 
of an electrically chargeco, spherically symmetric body. Discovered very soon [TJ [2] after the Schwarzschild 
black hole, the RN solution enjoyed constant interest since then, being studied from different directions. 
This includes solutions in the presence of a cosmological constant [3], and, generalizations to spacetime 
dimensions higher than d = 4 [4] . 

The stability of the RN black hole is an important issue. As discussed in [5], [5], the d = 4 solution is 
stable with respect to linear perturbations; the stability of higher dimensional solutions has been considered 
in [7], [S]. Moreover, a RN black hole is thermally stable for large enough values of the electric charge 

There is however a different type of instability that the RN solution can exhibit. This occurs when the 
Einstein-Maxwell (EM) system is considered as part of a larger theory, typically containing scalar fields. 
Then, the RN black hole can become unstable to forming hair at low temperatures. This results in new 
branches of hairy solutions of the larger theory, with the extra matter fields in a non- vacuum state. Usually, 
the hairy black holes are thermodynamically favoured over the RN solution. 

Although this type of instability was originally discovered for rf = 4 asymptotically flat, magnetically 
charged, RN black holes embedded in Einstein- Yang-Mills-Higgs theory [10], it has recently received a 
considerable attention in the case of black holes with anti-de Sitter (AdS) asymptotics. The main models 
considered in this context are the gravitating Abelian Higgs [TT] and Einstein- Yang-Mills [T2], [13] theories, 
augmented with a negative cosmological term. The resulting sets of hairy black holes possess Ricci-flat 



'^In d = 4 spacetime dimensions, due to the electric-magnetic duality, one may consider as well a magnetic charge source for 
this solution. 
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horizons, and have found interesting apphcations in the context of gravity/gauge duahty, providing a dual 
gravitational description of some superconductors (see e.g. [T3], [TS]; for a review of this domain, see [T^). 

The main purpose of this work is to discuss a a specific type of instability of the electrically charged RN 
black holes in d = 5 spacetime dimensions in context of the Einstein- Yang-Mills-Chern-Simons (EYMCS) 
theory. The presence of a cosmological term, positive or negative, is immaterial in the sense that we 
demonstrate qualitatively similar features in all these cases. This unified treatment of the various asymptotics 
is part of our brief here. 

Usually, the RN solution is written in Scwarzschild-like coordinates, with a line element 

ds^ = — — + r'^dnl - N{r)df, (1.1) 



N{r) 



with 



iV(r) = l--;^+e-, ™(0=A/o-^, (1.2) 

(where A = ^ejl? is the cosmological constant, with e = —1, 0, +1 for AdS, Minkowski and dS asymptotics, 
respectively, and d^\\ is the line element of the three dimensional sphere), and an C/(l) electric potential 

y(j)-y^-%- (1.3) 

Mq and Q corresponds (up to some constant factors) to the mass and electric charges of the solution; also 
we note that o? = 167rG/(3e^), with G the Newton constant and e is the non-Abelian (nA) gauge coupling 
constant. 

The instability discussed in this work occurs when embedding the RN solution ()l.ip - (|1.3|) in EYMCS 
theory, the Chern-Simons (CS) term in which plays a crucial role. In the absence of a cosmological constant, 
this model has been discussed in [T7], [TH], where the corresponding zero mode of the RN black hole was 
exhibited, together with the resulting set of hairy black hole solutions. The results there indicate the 
existence of a phase transition between RN solutions and the non-Abelian black holes, which generically are 
thermodynamically preferred. Solutions of the EYMCS model with a negative cosmological constant were 
discussed in [19j . 

One objective here is to argue that the instability of the RN black hole discovered in [ITl [18], [l9| 
for Minkowskian and AdS asymptotics respectively, occurs also in the presence of a positive cosmological 
constant in the action. As a result, a unified picture of this instability manifests itself in these three cases. 
For a fixed value of the CS coupling constant, a classical instability is developed if the horizon is sufficiently 
small. We find that a number of basic features of the resulting set of hairy black holes are the same both 
for asymptotically fiat, and, for (anti)-de Sitter solutions. Some of the EYMCS solutions with A < in [19] 
are interpreted here in this context. Finite mass solutions with dS asymptotics are also constructed. In all 
cases, an Abelian gauge symmetry is spontaneously broken near a black hole horizon with the appearance 
of a condensate of non-Abelian (nA) gauge fields there. 

The paper is organized as follows. In the next Section we define the model and propose a spherically 
symmetric Ansatz. Section 3 is devoted to the stability analysis of the RN solution in EYMCS theory. We 
derive the pulsation equation governing the evolution of radial nA perturbations and solve it numerically. 
Special emphasis is placed on the issue of marginally stable modes, an analytic expression being found in 
this case. In Section 4 we present numerical results and discuss qualitative properties of the hairy black hole 
solutions. For solutions with a positive cosmological constant, we construct new sets of black holes with nA 
hair approaching asymptotically the dS background. We conclude with Section 5, where the significance of, 
and further consequences arising from, the solutions we have constructed are discussed. 
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2 The model 



2.1 Field equations 

We consider the following EYMCS action 



1 1 



5=/ d^xV^{^^iR-2A)--Tr{T^,T^''}] - I d^'x Ccs, (2.4) 



M yl6nG 4 / JM 

(with J^^i, — dfj^Ai, — dyAp + e[ylp, Au] the gauge field strength tensor), and 



£CS = KenyparT^l- { A^ 



j:pvj:pa _ ^j:l^^j^Pj^o ^_ A"" A" A" 
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(2.5) 



is the CS term. The field equations are found by taking the variation of (|2.4p with respect to , Ap and 
are given e.g. in [19] . 

There is some arbitrariness in the choice of the gauge group, the only restriction being that it should be 
large enough to accomodate a nonvanishing electric potential, such that the RN black hole is a solution of the 
model. The electric potential is supported by the CS term, which would likewise vanish if the gauge group 
is not large enough. Subject to these criteria, the smallest simple gauge group supporting a nonvanishing 
CS terni is 50(6). 

In this work, as in [IF, we shall restrict to a 5*0(4) x 50(2) truncation of the 50(6) gauge group, i.e. 

A^ A-dr^adVY:, (2.6) 

where the magnetic potential A takes values in the 50(4) ~ SU{2) x SU (2) subgroup, while T is the Casmir 
operator corresponding to the generator of the 50(2) subgroup. 

For such a truncation of the gauge group, (|2.6p . it is important to note that the EYMCS model admits 
an equivalent formulation. Restricting to an SU(2) subgroup of 50(4), and discarding a boundary term 
arising in the CS term, one finds that (|2.4p reduces essentially to a Einstein- Yang-Mills-Maxwell model, 
with a Chern-Simons-type coupling term between the U{1) and nA fields 

S = J^ d^x^g {j^{R - 2A) - \fI,F^-' - \f^,r - -^e'''''«^^a^Fj^Fi,) , (2.7) 

where is the SU{2) non-Abelian gauge field (with / = 1,2,3 and field strength = i9^A^ — di^Aj^ + 
ssukA'^A^), and is the U{1) gauge field {f^i, = d^aiy — di,a^j, being the corresponding field strength). 
Variation of (|2.7p with respect to g^^, Aj^ and leads to the following equations: 

1 

r 



Rpn — -T^Rgpv + ^gpv = SnGT^^i,, 



d 

dxP- 
d 

^x^^ 

where 



Tpu = Fl^Fl^g'^P - -^g^^FipF^P' + fpa.Ug'^^ - ^gpuU^r^ (2.9) 

is the energy- momentum tensor. Some features of this model with A = have been discussed in |24) . where 
it was found that the value k = a/8 (with = 167rG/(3e^)) is special, being given by the coupling of 
the super- YM theory to the d — 5 supergravity. This allows to promote all d = 4 flat space instantons 
to extremal black hole and soliton solutions of the Eqs. (|2.8I) with a vanishing cosmological constant, see 
the discussion in [TH]. The model (|2.7p can also be thought as a truncation with a vanishing dilaton of the 
A/" = 4, d = 5 Romans' gauged supergravity model pS] . 

■^Then, for A < this model can be viewed as a truncation of the five-dimensional maximally gauged supergravity |20) . |21j . 
which plays a prominent role in the conjectured AdS/CFT correspondence |22| . |23l . One expects that the basic properties of 
the solutions in this truncation, persist in the full model. 
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2.2 A spherically symmetric ansatz and the Abelian solutions 

In this work we are interested in spherically symmetric solutions only, with a line element 

ds^ = -Na^dt^ '^1^^ r^dnl, (2.10) 

where N, a are functions of r and t in general. The gauge fields ansatz contains only two functions, a 
magnetic SU{2) gauge potential w, and an C/(l) electric one, V: 

= -(l + w(r,t))0^, a^-V(r,t)dt, (2.11) 

e e 

where 6^ are the invariant forms on S^, such that d9^ — ej^O'^ A 9'^ and dil^ = 6^9^ . 

The RN Abelian solution is recovered for N given by (II. 2p . a ^ 1, w = ±1, and V given by (II. 3p . This 
configuration has a black hole event horizon, located at r = r/j (which is the largest root of the equation 
N{r) — satisfying the condition iV'(r;j) > 0). The conserved quantities associated with the RN solution 
are the mass Ai = and the electric charge Q = ^f-Q. Other quantities of interest are the chemical 

potential $ = the Hawking temperature Th — jnr^C^ ~ "J^ + 2e-p-), and the entropy S = ^ (with 
Ah = 27r^r^ the black hole event horizon area). 

For the discussion of the solutions, it is useful to introduce the dimensionless quantities 

Q_ (2.^)V3Q _ rl A^ 

^ - A]i^ ' ^'''^ ^ " - (2,r2)2/3L2- (2.12) 

These parameters characterize uniquelly the RN black hole (asymptotically flat solutions have Y = {)). For 
example, several dimensionless quantities of interest in the thermodynamical description of the system can 
be expressed in terms oiU,Y: 



where /, tn and an are the reduced free energy, temperature and area, respectively. 



3 The instability 

3.1 Time dependent perturbations 

We now investigate the stability of the RN solution, when embedded in the EYMCS model. The functions 
iV, cr, w and V satisfy a complicated set of partial differential equations which are found after replacing the 
ansatz ()2.10p . ()2.1ip in the field equations ()2.8p . However, for the purpose of this paper, it is sufficient to 
study the equation for the magnetic gauge potential only : 

d ( dw\^ d f r dw\^ 2a-w{l — w'^) ^ ^ 2\^^ 



Also, in what follows, we use the scaling symmetries of the considered EYMCS system [T7], [15] to set a = 1. 

The instability is found for values of the nA magnetic gauge potential close to the vacuum value every- 
where. To consider fluctuations around the RN black hole, we only need to keep linear terms in 6w = w + 1 
(similar results are found for Sw = w — 1), SV = V — and SN = N ~ N^^-^\ Sa = a — 1. It turns out 

that the coupled equations separate and the equations for the metric components contains neither Sw nor 
SV. The instablity occurs if there are solutions of the form Sw = VF(r)e^*, with real fi. From p.l4p . one 
can see that W{r) is a solution of the equation 

+ N)W -n^W ^0, (3.15) 

dp^ dp 
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Figure 1. The square of the frequency Q,^ is plotted as a function of the Chern-Simons couphng constant k (left) 
and of the dimensionless parameter U = Q/rj^ (right). 
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Figure 2. The 'charge to radius' ratio U{k) — Qlr\ which characterizes the unstable Reissner-Nordstrom solu- 
tion where a branch of non-Abelian configurations emerges is shown as a function of the CS coupling constant k for 
solutions with a cosmological constant A < (left) and A > (right). For de Sitter solutions, we show also the 
reduced free energy / = (A^ — ThS)/Q of the critical Reissner-Nordstrom configurations as a function of reduced 
temperature tn ~ Th^/Q, the CS coupling constant there being the intrinsic parameter along that curve. 

where we have introduced a new 'tortoise' coordinate p defined by ^ = rN , such that the horizon appears 
at p — > — oo. Close to the horizon, the physical solutions of p.lSp behave as 

Wir) = woe^P, (3.16) 

with Wo a free parameter. As r — oo, the approximate form of the solution is 



w{r) — e 



(3.17) 



k=l 



for A = (where all coefficients Wk are expressed in terms of Wi, e.g. W2 = Wi(15 + 8Afoil^)/(8ri)), and 

oo 

^M=E5' (3-18) 



wi 



k=2 
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for A 7^ (again in terms of a single free coefficient W2)- 

The perturbation W{r) is determined by solving numerically the equation (j3.f 51) with the boundary 
conditions (|3.f 6l) - (|3.f 81) . by using a standard Runge-Kutta ordinary differential equation solver. In practice 
we use a shooting aproach in terms of the parameter fl^. 

In Figure I we plot the frequency as a function of the CS coupling constant k for several values of 
U — Q/rf^ (left) and as a function of U for several values of k (right). One notices that, given Q, r^, the RN 
black hole becomes unstable for all values of k > Kmim with — > as the minimal value of k is approached. 
Equivalently, for fixed k, a classical instability develops if the horizon is sufficiently small {i.e. ii U = Q/r\ 
is large enough). Similar results were found in the presence of a cosmological term. 



3.2 The marginally stable mode 

For = 0, Eq. p.lSp takes the simple form 

r(rW)'-4(l-^)Ty = . (3.19) 

The solutions of p.l9p are of particular interest since they indicate the emergence of a branch of static, 
hairy black holes with a nonvanishing magnetic field outside the horizon, w(r) ^ ±1. The existence of these 
configurations can be understood as follows [T7], [H]: the second term in p.lSp shows the presence of an 
effective mass term /i^ for W near the horizon, with ^ 1 — 8U{k). All physical solutions have /i^ < 0, 
with Q/rf^ — U{k) being a monotonic function of the CS coupling constant k. 

Although the equation (j3.19p does not appear to be solvable in terms of known functions, one can 
construct an approximate solution near the horizon and at infinity. As r — > r/j, the function W{r) behaves 
as 

= b + 2fcr.(r^-8.Q)L^ _ ^^^^ ^ _ ^^^^^^ 



2srl + {rt-Q^)L^) 
For large r, the approximate form of W{r) is 



W{r) ^^ + ^{AIo- AkQ)^ + . . . , for A ^ 0, (3.21) 
3e r° 



and 



7 2 7 1 
W{r) - - + —{Mo - it^Q)- + ..., for A = 0. (3.22) 

In the above relations, b and J are free parameters (in practice one sets b — 1 without any loss of generality)- 
Solutions interpolating between the asymptotics p.20|) and (|3.21|) . p.22|) are constructed numericalljij. 

In practice, the input parameters are k and Y = r\/Lp' . The "charge to radius" ratio U = Q/r\ results from 

the numerical output. In Figure 2 we show the profile of U{n) for several values of the cosmological length 

scale L (the event horizon radius is taken r/j = 1). 

An analytic expression of U{k) can be found by matching at some intermediate poini|f| the expansion of 

W{r) (and its first derivative) at the horizon, p.20p . to that at infinity, p.2ip . (|3.22p . Surprisingly, it turns 

out that the resulting expression 

[/ = = -4k + ^2 + 16k2 + 2eY, (3.23) 

provides a good approximation of the numerical data (typically with several percent error) and gives also an 
analytic explanation for some numerical results found in |17) for A = 0. For example, for large k, one finds 

^In this work we restrict our study to solutions with a monotonic behaviour of W{r). 
*In this approach, a new radial coordinate is introduced, x = r^/(l — x), with < a; < 1. 
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[/ ~ (1 + eY)/{4K), with Uk =1/4 for A = 0. Also, the condition for tachionic mass gives a minimal value 
of the coupHng constant k 

^ ^ 8Vl + 2ey' ^^'^^^ 

which imphes k > 1/8 for A = 0. 

The knowledge of the function U allows us to reconstruct, after replacing in (I2.13|) . the thermodynamical 
parameters of the unstable RN solution as a function of k,Y. For example, when embedded in EYMCS 
theory with A = 0, the RN black holes possess a marginaly stable mode for a value of the dimensionless ratio 

TrAh = ^ (l - (4« - ^2 + 16^2)2) ' , (3.25) 

while tn = ThVQ ^ as k ^ 1/8. 



4 Black holes with non-Abelian hair 

The instability of the RN solution pointed out in the previous section can be viewed as an indication of 
the existence of a branch of non-Abelian solutions with nontrivial magnetic non-Abelian fields outside the 
horizon. These configurations appear to be the natural end points to which the instability of the RN solution 
leads. 

Such solutions are constructed nonperturbatively by solving numerically the set of equations following 
from (121]) : 

2aw{w^ - 1) 2 IN f ^^^' v / 2 ^^ , 

iraNw) = |-8kI/(w —1), I I = 24k(w — Ijw , 

r a 

where we set = 167rG'/(3e2) = 1 and define m(r) according to 

The function m(r) is (up to a constant factor) the local mass-energy density. 

With the scale fixing adopted, the reduced equations depend on two input parameters ~ the (A)dS length 
scale L and the CS coupling constant k. The invariance of (I4.26P to changing the signs of w, V and k allows 
to study solutions with strictly positive values of k only. 

We notice that the equation for Vir) has the first integral [18] 

T^' = ^(Q + 4/t(w-2)(w + l)2), (4.28) 

with Q an integration constant which is taken to be positive without any loss of generality. 

For any sign of A, the expression of the solutions as r — oo as imposed by the finite mass/energy 
requirement, i^j 

m(r) = M-^-t^J^ + ..., a(r) = l-^ + ..., «;(r) = -1 + ^ + . . . , F(r) = - ^ -I- . . . , (4. 
where J, M, Vb are parameters given by numerics which fix all higher order terms. 



'Choosing ui(oo) = 1 leads to similar results. 
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We want the line element (|2.10p to describe a black hole, with an event horizon located at r — > 0- 
In the vicinity of the event horizon, one can write a formal power-series expansion of the solution 



m{r) Too + ^ mk{r - , cr{r) = cr/i + ^ crfc(r - rh)'', 

k=l k=l 

oo oo 

w{r) =Wh + ^Wk{r-rh)'', V{r) ^J^Mr - rn)'' ■ (4.30) 
fe=i fe=i 

All coefficients in the above series are determined through reccurence relations by Wh and ah- One finds e.g. 
vi = '^iQ + 4k(1 + Wh)^iwh ~ 2)), 

'"i = TTs + 32kQ(1 + Whf{wh - 2) + ^^n^iwn - 2)^(1 + Whf + 3r^(M;^ - 1)^), (4.31) 



2l?{^KrhV\ + crfeW/i)(wf^ - 1) 



Thus this expansion (assuming that its radius of convergence is nonzero) defines a two-parameter family of 
local solutions labeled by the values of w{r) and <j{r) on the horizon. 

The only conserved quantities associated with these solutions are the mass M. and the electric charge Q 

37r 47r^ 

M^^M, Q=—Q. (4.32) 
8G e 

Other quantities of interest are the chemical potential $ = — , the Hawking temperature Th , and the entropy 
S 

T,^l^airh)N'irh), ^ = ^ = (4.33) 

One should note that, similar to other uA hairy black holes, there is no conserved quantity associated with 
the parameter J which appears in the large-r asymptotics (|4.29p . 

The dS solutions have some special features due to the presence of a cosmological horizon at r ~ rc > ■ 
(See the discussion in Section 4.) 



4.1 The cases A < 

The asymptotically flat solutions of this model have been discussed at length in [TS]. The results there 
indicated that the nA solutions emerge as perturbations of the RN black holes for a critical charge to radius 
ratio U, which is fixed by the CS coupling constant k. Hairy black holes appear to exist for any value of 
K > 1/8. For a fixed r^, the value of the electric charge parameter Q is also important, some basic features 
of the solutions depending on whether Q is less or greater than 

Q(") = 16k. (4.34) 

For fixed k, a plot of the reduced horizon area an vs. the reduced Hawking temperature tn reveals the 
existence of several branches of solutions. For Q ^ the picture is rather similar to that valid for RN 

black holes [18] . The temperature reaches its maximum at some intermediate value of the event horizon 
radius, an extremal black hole being approached for a minimal value of > 0. The lower branch of solutions 
has positive specific heat. Solutions with Q = Q^'^'> are special, since one finds a single branch of thermally 
unstable configurations. These solutions emerge again from a critical RN solution with rh = \/ Q^'^^ /U{n) 
and can be continued for an arbitrarily small value of the event horizon radius. As rh — > 0, these black holes 
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Figure 3. The reduced event horizon area an = Ah/Q^^'^ (left) and free energy / = [M ~ ThS)/Q (right) 
are plotted vs. the reduced temperature tu = ThQ^^^ for AdS branches of non-Abelian solutions with several values 
of the electric charge and a given value of k. 



3 



2.5 
2 

a 

^ 1.5 

1 

0.5 



"0.6 0.75 0.9 1.05 ^ 1.2 1.35 1.5 

Figure 4. The order parameter J which enters the asymptotics of the magnetic gauge potential at infinity is 
shown as a function of the Hawking temperature Th for a fixed value of the Chern-Simons coupling constant k and 
several values of the electric charge. Note that Th is normalized with respect the temperature of the critical Abelian 
solution where a branch of non-Abelian solutions emerges as a perturbation. 

approach a set of globally regular particle-like solutions. More details on the properties of the A = 
configurations togheter with typical profiles can be found in |17) . |18) . 

The basic features of the AdS solutions were discussed in [19]. However, the main emphasis there was on 
the solutions of the full 5*0(6) model, featuring two extra gauge potentials apart from w, V. The SU{2) x U (1) 
solutions of interest here corresponds to the family of solutions which was referred in [12] as the 'fundamental 
branch of EYMCS solutions'. We have reexamined the properties and the basic features of these solutions 
and have found that they emerge indeed as perturbations of the RNAdS configurations. The solutions exist 
for values of k, above some minimal value which is well approximated by p.24p . 

We did not investigate in a systematic way the thermodynamics of the solutions corresponding to A 7^ 0. 
(Recall that the presence of a negative cosmological constant leads to a complicated thermodynamical picture 
already for RN black holes [2S], [13 )■ The set of solutions investigated so far have a relatively small ratio 
Y = rf^/L"^ and reveal a picture which is rather similar to the asymptotically flat case. For Q 7^ Q*-^-*, a 
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temperature-area plot indicates the existence of two branches of solutions (see Figure 3 for typical results). 
The upper branch of solutions emerges from a critical RN black hole with — ^ Q IU{k) and is thermally 
unstable, i.e. the temperature increases for decreasing r^. The lower branch always possesses a positive 
specific heat, the Hawking temperature vanishing there for a minimal value rj^™'"^ of the event horizon 
radius. As r — >■ r|j"*"'*, a zero temperature non-Abelian black hole solution with a regular horizon is 
approachecfl. As noticed already in [12], the solutions with Q = Q'-'^^ = 16k exist for arbitrarily small values 
of r/i, a particle-like solution with a regular origin being approached as rh — > 0. 

In Figure 4 we plot the dimensionless quantity J /Q (with J the order parameter which enters the first 
relevant term in the large— r expansion of the magnetic gauge potential) as a function of the reduced Hawking 
temperature and several values of the electric charge parameter. One can see that for Q < Q^'^\ J approaches 
a constant, nonvanishing value for small enough values of T^. 

In the discussion of both asymptotically fiat and AdS solutions, we have restricted to configurations 
with a monotonic behaviour of the magnetic gauge potential w(r). Solutions with local extrema of w{r) do 
exist, but it is likely that they are always thermodynamically disfavoured because spatial oscillations in w 
increase the total mass. Also, as expected (since they start as perturbations around RN black holes), for 
both Minkowski and AdS asymptotics, we have found nA solutions where w{r) has no nodes. 

4.2 A > 0: de Sitter solutions 

Similar solutions should exist for positive values of A, as indicated by the presence of an unstable mode also 
in that case. However, the situation is much more involved in this case due to the presence of a cosmological 
horizon. 

The dS configurations are found again by solving numerically the Eqs. (j4.26p . Some of the basic features 
of these solutions are generic for any black hole in a dS background. For A > 0, the solutions possess a 
cosmological horizon located at r = Vc > rh. In this case, one needs another set of boundary conditions at 
the cosmological horizon, with N{rc) = 0, N'{rc) < 0, a{rc) > 0. A formal power series similar to (|4.30|) 
can also be written at r = Vc, the only free parameters being w{rc) and cr(rc). Moreover, both the event 
and the cosmological horizons have their own surface gravity (and hence different temperatures.) Outside 
the cosmological horizon r and t changes character {i.e. r becomes a timelike coordinate for r > Tc). A 
nonsingular extension across this null surface can be found just as at the event horizon of a black hole. The 
regularity assumption implies that all curvature invariants at r = rc are finite. Also, all matter functions 
and their first derivatives extend smoothly through the cosmological horizon, e.g. in a similar way as the 
U{1) electric potential of a RNdS solution. 

The mass and electric charge of the solutions are computed outside the cosmological horizon, as r — >■ cx) 
(note that the asymptotic expansion (14.29^ still holds for A > 0, although N{r) < for r > Tc). Similar 
e.g. to the Schwarzschild or RN black holes, the total mass, as found by employing the quasilocal formalism 
in [31], is = —|^M (thus with a minus sign as compared to the definition (|4.32p . which is well-known 
feature of the dS solutions). 

The field equations were solved in this case by employing a collocation method for boundary-value 
ordinary differential equations, equipped with an adaptive mesh selection procedure [32) . As in the Einstein- 
Maxwell case, the cosmological horizon radius rc is a function of A. In practice, we solved the equations first 
on the interval [rh,rc\, choosing Vc by hand and imposing regularity conditions of the solutions at r = rc. 
This allows to determine the values of the functions (m, a, w, w' ,V) at r — rc, as well as the numerical value 
of A corresponding to the choice of cosmological horizon. The integration on [rc, oo) can then be performed 
as a second step, leading to the knowledge of the solution on the full r— interval. 

Our numerical results indicate the existence of nontrivial solutions of the equations (14.26^ approaching 
as r — > cxD the dS background. The profile of a typical black hole solution with nA hair is shown in Figure 5. 
One notices that the Yang-Mills fields do not vanish at r — and persist beyond the cosmological horizon. 
Thus these solutions may be thought of as cosmic colored black holes. Also, similar to the A < case, 
solutions with a regular origin are found for the value of the electric charge Q — IQk only. 

^ The existence of these extremal configurations is supported by an exact AdS2 x solution of the Eqs. 114.261 1 ^19^. This 
exact solution describes the near horizon geometry of an extremal black hole with nA hair. 
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-0.5 0.5 1 1.5 



Figure 5. The profile of a typical hairy black hole solution with de Sitter asymptotics is presented as a func- 
tion of the radial coordinate r. m{r) and a{r) are metric functions, while V{r) and w{r) are electric and inaguotic 
gauge potentials, respectively. 
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Figure 6. Several relevant parameters are plotted as a function of the magnetic gauge potential on the horizon 
(left) and of the electric charge parameter Q (right) for de Sitter black hole solutions with a fixed value of the 
Chern-Simons coupling constant k and a given event horizon radius. 



However, determining the pattern of the dS solutions in terms of the input parameters k, L represents 
a very complex task which is outside the scope of this paper. Instead, we have studied particular sets of 
solutions for several choices of k, L only. Several parameters characterizing these solutions are presented on 
Figure 6. There we show first the dependence of the solutions on the value of the magnetic gauge potential 
on the horizon for fixed values of the CS coupling constant k, event horizon radius r/j and cosmological length 
scale L. The quantities plotted there are normalized with respect to those of the critical RN-dS solution 
which possess a marginally stable mode. The dependence of the solutions on the value of the electric charge 
is also shown in Figure 6 (right) for a different set of relevant parameters. 

When studying the pattern of solutions for a fixed value of the electric charge Q and increasing the event 
horizon radius , we find that generically the branch of nA black holes ends up in an extremal solution for 
r/i — r^. However, the picture is quite complicated and seems to depend in a nontrivial way on the value 
of K. For example, for some values of k, wc have found that when Wh becomes sufficiently small, the metric 
function N{r) develops a local minimum at some intermediate value of the radial variable, say r = ro > rh 
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and tq < re- The value -/V(ro) tends to zero for some minimal value Wh = Wh^c and the solution therefore 
bifurcates into an extremal RN black hole in the limit Wh — Wh.c- Inspecting the matter fields, we further 
note that, for the limiting solution, the gauge field function i/;(r) is not constant for r g [f'/i,fo)- However, 
the magnetic field vanishes for r £ [ro,oo), i.e. uP'^r) = 1, and at the same time, the electric potential V 
approaches the one corresponding to a RNdS black hole. As a consequence, an extremal RNdS black hole 
is formed for r G [ro , oo\ . This kind of bifurcation resembles a phenomenon observed some years ago in the 
context of gravitating monopole and black holes [33]. The control parameter is there a — Gjv^ (with G the 
Newton's constant, and v the expectation value of the Higgs field). It was shown that, when a becomes large 
enough, the metric function g„ (r) also develops a local minimum at some intermediate value r = rg and the 
gravitating monopole (or monopole black hole) bifurcates into an extremal RNdS solution for re [rg , cxd] . 
Here the phenomenon occurs in the absence of a Higgs field, the matter field is purely Yang-Mills type and 
the control parameter is k, which like the Higgs VEV, is a dimensionful quantity. 

We conjecture that, similar to the AdS case, a positive cosmological case allows for the existence of a 
larger set of spherically symmetric solutions, with four non-trivial gauge functions (i.e. for the full group 
S'0(6)). As argued in [1^, these solutions are unlikely to exist in the asymptotically flat case. 

We close the discussion here by noticing that for k, = 1/8, by using the generation technique proposed in 
[25] . one finds the following exact solution 

2t/L 7 _ 

= + '-'^"a) - ^'(^, Orfi', <r) = V{r, t) = F{r, t), (4.35) 

where F{r,t)~^ = 1 + e^^*/-^ + (r^i~77^) ^^'^ ^ ^ This solution is written in an inflationary 

coordinate system and describes a non-Abelian deformation of the d = 5 extremal RN-de Sitter black holes 
discussed e.g. in ['29 . The time-dependence enters in a natural manner in this solution and for large r, the 
dS metric in planar coordinates is approached. To some extent, the properties of the nA solution (j4.35l) can 
be discussed in a manner similar to the case of RNdS black hole written in planar coordinates, see e.g. [29j . 
In this approach, the mass is computed at far past or far future infinity (depending on the sign of L) by 
using the quasilocal tensor of Brown and York [30] (augmented by suitable boundary counterterms [31]), the 
resulting expression being A4 = — Stt^. 



5 Further remarks and conclusions 

In this paper we have presented arguments that the d — 5 electrically charged RN black holes become 
unstable when SU{2) nA fields are introduced in the action, with a specific coupling to the U{1) field. This 
instability allows for dressing RN solutions by nontrivial configurations of nA magnetic fields that arc not 
associated with a Gauss law, thereby leading to an evasion of the no-hair conjecture. Although some of the 
resulting nA solutions, with A < and A = 0, were already displayed in the literature [T7], [H], we have 
presented here a unified framework, valid for both asymptotically flat and (anti)-de Sitter solutions. In this 
context, new nA solutions approaching a dS background at infinity are presented and discussed. 

The solutions constructed here all have finite energy by virtue of the presence of the Chern-Simons 
(CS) term, whose scaling properties are adequate for that purpose. In this sense, employing a CS term 
is an alternative to the use of higher order Yang-Mills curvature terms, irrespective of the value of the 
cosmological constant [34j . For the dS case, a higher order YM curvature term was used in [351 to obtain 
finite energy nA solutions. The use of a CS term is more restrictive than that of higher order YM curvatures, 
since in any given dimension there is a unique (single trace) CS density. By contrast, the solutions in that 
case present a richer pattern. 

Our study reveals that some of the features of the solutions are rather generic and independent of the 
presence of a cosmological term in the action. In all cases, the nA black hole solutions emerge as perturbations 
of the RN solution, which becomes unstable when embedded in a larger gauge group. Moreover, for any 
asymptotics, a set of particle like solutions is found for a critical value of the electric charge, namely onlj[3 
for Q = 16k. 

This feature can be understood by noticing that the cosmological constant does not enter the first integral 1 14.281 1. Since 
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Following the approach in [T5] , we could prove that some of our AdS solutions are also stable with respect 
to first order of perturbations. The equations for A < are rather similar to those given in Section 4 of 
[15] for flat asymptotics, so we have not display them here. As expected, all stable solutions we have found 
have a nodeless magnetic potential w{r). Physically, this means that the ballance between the attractive 
gravitational force and the repulsive gauge U{1) field is not affected when adding small enough non-Abelian 
fields. The corresponding question of stability of dS solutions is complicated by the presence of a cosmological 
horizon and we defer it to future work. 

It may be interesting to contrast the situation with the AdSs solutions of the pure EYM system {k = 0) 
discussed in [36l [37]. In both cases, we regard a U{1) subgroup of the YM fields as the gauge group of 
electromagnetics, and persuade the off diagonal gauge bosons to condense outside the horizon. However, 
in the absence of a CS term, EYM black holes with finite mass exist only in AdS spacetimes, which are 
characterised by a Ricci flat horizon. Besides, such black hole solutions do not admit particle-like, globally 
regular limits. In contrast, the EYMCS model discussed here admits also finite mass AdS solutions with a 
Ricci fiat horizon [38] . in addition to the solutions with a spherical horizon. 

The EYM-A black holes in [36] describe, via gauge-gravity duality, a p— wave superfluid state in a strongly- 
coupled conformal fleld [37]. The possible relevance of the EYM solutions with a CS term in providing 
useful analogies to phenomena observed in condensed matter physics remains an interesting problem to be 
understood in the future. At least in the AdS case, we expect our solutions to provide duals of spherical 
superconductors. Further progress in this direction requires a computation of the conductivity as a function 
of frequency by perturbing the YM fields, a problem which is currently under investigation. It should be 
stressed that our results here signal the potentially important fact that some basic features exhibited by 
such solutions do not require the presence of a cosmological term in the action. This may suggest that the 
solutions with fiat or dS asymptotics may also be useful in describing certain aspects of condensed matter 
systems. 
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